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The phenomenological analysis of fully spin-polarized quantum Hall systems, based on holomor- 
phic modular symmetries of the renormalization group (RG) flow, is generalized to more complicated 
situations where the spin or other "flavors" of charge carriers are relevant, and where the symmetry 
is different. We make the simplest possible ansatz for a family of RG potentials that can interpolate 
between these symmetries. It is parametrized by a single number a and we show that this suffices 
to account for almost all scaling data obtained to date. The potential is always symmetric under 
the main congruence group at level two, and when a takes certain values this symmetry is enhanced 
to one of the maximal subgroups of the modular group. We compute the covariant RG /3-function, 
which is a holomorphic vector field derived from the potential, and compare the geometry of this 
gradient flow with available temperature driven scaling data. The value of a is determined from ex- 
periment by finding the location of a quantum critical point, i.e., an unstable zero of the /3-function 
given by a saddle point of the RG potential. The data are consistent with a G K, which together 
with the symmetry leads to a generalized semi-circle law. 
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I. INTRODUCTION 

The low energy physics that emerges at large scales from 
strongly interacting electrons in a background magnetic 
field and confined to two dirty dimensions, i.e., the quan- 
tum Hall effect (QHE), is not accessible by the conven- 
tional perturbative expansion of the local microscopic 
theory. Similar problems arise in high energy physics, 
where self-interacting gauge bosons give rise to an in- 
soluble set of highly non-linear differential equations. 
One of the earliest and most successful methods de- 
veloped to circumvent this problem used to be called 
"phenomenological" -, a precursor of what today is called 
effective field theory (which includes all quantum field 
theories in need of an ultraviolet cutoff). The essence of 
this idea is to use some of the global properties of the 
theory — in the case of chiral models (i.e., low energy 
QCD and current algebra) the geometry of the global 
symmetries that are observed in the hadronic spectrum. 
These are always approximate and may be broken in 
some phases, giving rise to pseudo-Goldstone bosons (the 
pions) rather than massless bosons. 

The "phenomenological" approach to the QHE is sim- 
ilar in spirit, since it exploits emergent symmetries to 
harness the phenomenology of the system, but the sym- 
metries are very different. Since these are found to be 
infinite, discrete, non-abelian and holomorphic, a new 
set of mathematical tools is required, as well as a pro- 
foundly geometric understanding of the renormalization 
group (RG) that we now describe. 



A. Renormalization 

An RG flow is a vector field on the space of those pa- 
rameters that are relevant at the chosen energy scale. 



The QHE is parametrized by the conductivity tensor 
CT*-' , or equivalently its inverse, the resistivity tensor 
p^^ = {(j~^y^ . These transport tensors are non-trivial be- 
cause the background magnetic field breaks parity (time- 
reversal) invariance, whence the off-diagonal Hall coeffi- 
cient is permitted by the generalized Onsager relation. 
A conventional Hall bar exhibiting fractional plateaux is 
usually a high mobility (/i '-^ IQ^ — lO^cm^/Vs), low den- 
sity (n ~ 10^ — 10^^/cm^) GaAs/AlGaAs heterostructure 
that confines all charge transport to a single isotropic 
layer of size x Ly, aligned with a current / in the x- 
direction so that R^,x — {L^, /Ly)p*^ . The Hall resistance 
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-p^^ is quantized in the fundamen- 



tal unit of resistance, /i/e = 25.812807557(18) kQ, while 
the dissipative resistance i?^ = p^^/D is rescaled by the 
aspect ratio □ = Ly/L^- 

We choose to label the low-energy parameter space by 
the complexified resistivity p = p^^ +ip^^ = —p^j +ipjj, 
or equivalently by the complexified conductivity a = 
a^y + ia^^ = a^j + «cr^. These complex coordinates 
take values in the upper half of the complex plane, 
H(o') = {cr G C|Scr = CT^ > 0}, because the dissipative 
conductivity (resistivity) is positive. The reason for not 
including the real line in our definition of the parameter 
space will soon become clear. 

The tangent vectors and to the RG flow on the 
two dimensional space of conductivities are the physical 
(Gellman-Low) /3-functions, which measure how fast the 
parameters (couplings) renormalize when the scale pa- 
rameter A changes: 



^ dt ' 



/3^ = ^. 

' dt ■ 



t = ln(A/Ao) 



Very little can be said in general about the proper- 
ties of RG flows, except that the topology of the flow is 
determined by the fixed points. The fiow ends at infra- 
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red (IR) stable fixed points (©), which in the quantum 
Hall case are the plateaux observed at rational values 
— € Q of the complexified conductivity. This set 
P C Q of IR fixed points are therefore the only real values 
that should be included in the physical parameter space, 
H0 = H U P. Furthermore, physical (contravariant) /3- 
functions have simple zeros at quantum critical points 
tT(g) € H0 for the localization-delocalization transition. 
Every member of this set E C of critical points must 
be a proper saddle point of the fiow, i.e., there should be 
both attractive and repulsive directions. 

In Ref. 2 it was proposed that a lot more can be said 
about RG flows in the QHE. Careful examination of the 
fixed point structure (i.e., the set PUE of stable and semi- 
stable fixed points) probed by scaling experiments reveal 
that there appear to be emergent symmetry groups T 
acting on the parameter space by fractional linear trans- 
formations 7(cr) = (acr -|- b)/{ca + d). The group ele- 
ments 7 € r are given by integers a, b, c and d satisfying 
det 7 = det(a, b;c,d) — ad ~ be — 1, and some additional 
constraints that distinguishes the different sub-groups of 
the full modular group PSL(2,Z). These so-called mod- 
ular symmetries can be used to constrain the RG flow, 
rendering them essentially unique in maximally symmet- 
ric cases. 

In the absence of a rigorous derivation of the emergent 
symmetry T from indisputable microphysics we have no 
other guide than experiment to aid us in the identifica- 
tion of r. The full modular group is the simplest con- 
ceivable emergent symmetry, but it is not of direct phys- 
ical interest because it does not agree with data, as we 
shall discuss at length below. It is, however, the start- 
ing point for the mathematical discussion of these sym- 
metries. Our presentation mirrors the conventional one 
found in mathematics, which systematically develops the 
theory of modular transformations from the group theory 
of PSL(2,Z) and its subgroups. This may seem unnatu- 
ral from a physics point of view, since PSL(2,Z) is not 
a physical symmetry, but appears to be the only way to 
obtain a clear picture of these symmetries and how they 
are related. This is crucial in our work, since we are 
looking for a framework that interpolates between the 
different effective theories that are relevant for systems 
with different low energy degrees of freedom. 



B. Modular symmetry 

This type of symmetry is not without precedent in con- 
densed matter physics. Kramers- Wannier duality is a 
discrete Z2-symmetry Dj^^y that acts on the parame- 
ter space of the classical Ising model. It swaps the 
high temperature phase with the low temperature phase, 
which means that it is its own inverse, Dkw = ^kw 
^KW ^ ^- This is the prototype of the duality transfor- 
mation that makes modular symmetries interesting. 
The first and most important observation^ is that any 



F-symmetry partitions the parameter space into univer- 
sality classes, with each phase "attached" to a unique 
(plateau) fixed point ® on the real line. This follows from 
the mathematical fact that compactifying the topology of 
the space on which a modular group acts gives precisely 
the physical parameter space H®. Hall quantization is 
therefore an automatie and unavoidable eonsequence of 
modular symmetry. This is the first example of a re- 
markable confluence of quantum Hall physics and mod- 
ular mathematics. 

These quantum symmetries only emerge at tempera- 
tures that are so low, compared to the scale of the effec- 
tive low energy modes of the quantum condensate, that 
thermal fluctuations are swamped by quantum fluctua- 
tions. So, unlike space-time or gauge symmetries, these 
global discrete parameter space symmetries are not exact, 
and we need to know how approximate they are. The ap- 
propriate theoretical framework for addressing this ques- 
tion is effective fleld theory (EFT), which is a systematic 
expansion in the inverse mass scale that separates the 
microscopic physics (here quantum electro-dynamics in 
a disordered medium in a strong background field) from 
the large scale physics being probed by transport exper- 
iments. It is presumably the leading term in this expan- 
sion that is F-invariant. Higher order terms may break 
the symmetry, but are suppressed by inverse powers of 
the mass gap. If F were exact the locations of all fixed 
points would be given by rational numbers, so the prox- 
imity of the experimental values to the predicted rational 
values is a very sensitive test of the symmetry. Two ex- 
periments of this type are shown in Fig.[T] 

Fig.[l]^a) shows that universality, which is the key con- 
cept powering the RG approach to EFT, has been verified 
in the QHE to an unprecedented accuracy^. Any non- 
universal contribution to the plateau value was found to 
be less than a few parts per trillion. This is the reason 
that the QHE will be used as one of the corner stones 
in the new SI system of metrology under construction. 
These experiments have also found that the plateaux val- 
ues are rational at the per billion level. 

If we interpret this as experimental evidence for a mod- 
ular symmetry, then other consequences of modular sym- 
metries should also be very accurate and provide further 
tests of the modular model. For example, the fixed point 
set of each F cannot be manipulated: all the plateaux 
and quantum critical points following from F must be 
includecP. As soon as one quantum critical point is 
pinned down the location of all the others is fixed by 
the symmetry. In maximally symmetric cases, which in- 
cludes the spin polarized system, there is no freedom at 
all: the exact location of all fixed points follow directly 
from the symmetry. Only recently have scaling experi- 
ments probed sufficiently low temperatures for this idea 
to be properly tested, even in the simplest case of the 
fully spin polarized QHeP. 

Fig.[l|^b) compares an experimeniP determining the lo- 
cation of the quantum critical point in the delocalization 
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Figure 1: (Color online) Two experimental tests of the univer- 
sal data encoded in a modular symmetry, (a) Measurements 
of Apjj = pj^(GaAs/GaAlAs) — p^(Graphene) for the second 
Hall plateau = = 1/2 [h/e^] has verified universality 
at the per trillion level (adapted from Ref.[3|). (b) The loca- 
tion of the 3® 4 quantum critical point, predicted by modular 
symmetry to be at p^ = {7 + i)/25 [h/e^] (the relevant part of 
the modular phase diagram is shown in the inset j^, has been 
verified experimentally at the per mille level (adapted from 
Ref.O. 

transition between the the third and fourth Hall plateaux 
(3 iSip 4 in our notation) , with the location predicted by 
the em ergent modular symmetry identified two decades 
agcPEl. With the lowest temperatures around 10 mK 
we expect that any modular symmetry is fully emerged 
and controlling all universal aspects of transport activity. 
The agreement at the per mille level inspires further con- 
fidence in the modular idea. Observe that the only trace 
that appears slightly displaced from the theoretical value 
is obtained at the highest temperature (~ 0.5 K), and 
that it exhibits a very poorly articulated Hall quantiza- 
tion. This suggests that the system is not deeply into 
the quantum domain when the temperature exceeds a 
few hundred milli-Kelvin, and therefore is in a regime 
where the modular symmetry is expected to fail. It is 
therefore rather surprising how accurate it still appears 
to be. Given a proper understanding of the separation 
of scales in this system, it would be possible to estimate 
how accurate the EFT is, and therefore to quantify how 



approximate the symmetry is as a function of scale. 



C. Modular RG flow 

Since the actions of the renormalization group and the 
modular group must be consiste nt, F-symmetry forces 
the RG flow into a straight-jacketP^. Not only are the 
RG fixed points, including the quantum critical (saddle) 
points, mapped into each other, so are the /^-functions. 
This means that the flow rates everywhere in one phase 
are mapped into the flow rates in any other phase, and 
the flows observed experimentally in the spin-polarized 
QHE appear to satisfy this prediction as well'*. In partic- 
ular, it follows that the critical exponents X±, which are 
the (inverse) flow rates in the principal directions near 
a critical point, are "super- universal" : they are always 
the same, independent of which quantum phase transi- 
tion is considered^. This is consistent with measure- 
ments of the relevant exponent (A+ > 0) for different 
Hall transitionJ^!^ as well as numerical "experiments"'^. 
The irrelevant exponent (A_ < 0) appears to be all but 
impossible to measure, but numerical work reveals the 
striking possibility that A_ ~ — A+~. 

This is evidence of an analytic structure that has been 
used to identify the exact geometry of the RG flowlSllo] 
We can combine the two real /3-functions into one com- 
plex function (3'^ = + if3'^, but this /3-function can 
not be holomorphic since holomorphic (complex analytic) 
functions do not have proper saddles, as required by 
physics. This is seen by expanding a holomorphic func- 
tion near a (simple) vanishing point zq. There is only 
one real eigenvalue A = A+ = A_, which cannot vanish 
since zq would then not be a simple zero. Consequently, 
Zq can be a source (A > 0) or a sink (A < 0) for the flow, 
but not a saddle point .1^ 

In another example of the convergence of quantum Hall 
physics and modular mathematics, this is not an addi- 
tional constraint, since holomorphic modular contravari- 
ant (physical) vector fields do not exist at all. This fol- 
lows from the observatiorP that P"' transforms as a con- 
travariant vector field when a ^ a' = 7(cr) under a mod- 
ular transformation 7 e F: 

/3'^'= (^)/3^ = (cc7 + d)-2/3-, 

where the modularity constraint ad — be = 1 has been 
used. A non-singular holomorphic function that trans- 
forms like this is called a modular form of weight w = —2, 
and the most basic fact of modular mathematics is that 
no form with negative weight exists. In other words, 
physics and mathematics are in agreement that such 
functions should not exist. 

The spectrum of holomorphic modular forms is very 
sensitive to the choice of symmetry group F. If we want 
the form to be covariant under the full modular group 
PSL(2,Z), then the lowest possible weight \s w — A. For 
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the sub-groups found to be of relevance to the QHE the 
lowest weight is w = 2. Since the covariant /3-function 
has this weight, 

Per can be both modular and holomorphicW The combi- 
nation of these constraints gives an extremely rigid struc- 
ture to the EFT, without which we would have no hope 
of proceeding along this path. Real modular functions 
exist, but are all but impossible to work with. 

If the EFT does have a complex structure respected by 
an emergent modular symmetry, then we are in unprece- 
dented circumstances that allow the exact determination 
of global properties of the RG flow. There are two inde- 
pendent arguments favoring this, one experimental and 
one theoretical. 



D. Holomorphic RG flow 

Observe first that a physical parameter space should be 
an ordinary Riemannian manifold with metric G, which 
in our context must also be F-symmetric, i.e., a real 
modular form of weight {w^w) — (2,2). The natural 
geometry of is hyperbolic and the canonical met- 
ric is the Poincare metric. This is the Kahler metric 
with components Gh = Gh = dad^K = l/cr%, where 
K ~ ln((7£)|/((T)p) is a F-invariant Kahler potential and 
/ is, up to a phase ("multiplier system"), a holomorphic 
modular form of appropriate weight (w — 1/2). K is a 
physically reasonable potential, in the sense described in 
this section, if / is the Dedekind 77-function, but we do 
not need detailed knowledge of the metric here. Since we 
are assuming that a non-singular EFT exists for finite 
values of a, this metric is non-singular at critical points 
and therefore invcrtible, and we have: 

^phys ^I3<^ ^ G'^'^p^ . (1) 

For finite values of a we can therefore quarantine non- 
holomorphicity of the physical /3-function to the metric. 
So-called "holomorphic anomalies" may appear at sin- 
gular values of cr, but they will not concern us here. 
If /3o- is holomorphic the pseudo-experimental fact that 
A_ — —\+ follows from eq. ([T]), no matter what non- 
singular value the physical metric takes at the critical 
points. 

A second argument in favor of a holomorphic /3o- follows 
from our expectatiorP that this vector field is a gradient 
fiow. This means that 

/?, = , (2) 

where d„ — d/da and the RG-potential $ G M is a kind of 
"vacuum entropy" that counts the number of (effectively) 
massless degrees of freedom at critical points. 



The existence of these potentials has been proven for 
two-dimensional unitary quantum field theorieJ^ where 
they are known as "C-functions". Quantum Hall dy- 
namics is essentially two-dimensional since the phase 
space of incompressible quantum fluids is effectively two- 
dimensional (the spatial coordinates are canonically con- 
jugate, as in Onsager's vortex dynamics). Furthermore, 
a similar result is expected to hold in any dimension, so 
we may reasonably expect that the theorem applies. 

This so-called "C-theorem" guarantees that an RG po- 
tential exists that smoothly interpolates between confor- 
mal fixed points, which are critical points of the poten- 
tial, and that the physical /3-function is obtained from 
the gradient vector field generated by this potential us- 
ing a metric on parameter space that will not concern us 
here. Suffice it say that it can be calculated directly from 
the EFT, if this is known. Since it is smooth and positive 
definite it does not affect the topology of the RG flow, 
only the absolute values of the flow rates. By construc- 
tion this /3-function is completely normal: it vanishes at 
critical points (where the tangents of the potential are 
flat), and the critical exponents arc given by the princi- 
pal curvatures of the potential at these saddle points. 

While Zamolodchikov's proof is quite explicit, compu- 
tation of his RG potential requires access to various corre- 
lations functions, i.e., essentially the full effective action. 
Since this is not available for the QHE a more oblique 
approach is needed. This is provided by Friedan's proof 
of the C-theorerrfS]^ which uses only general properties of 
spectral functions. This proof shows explicitly that the 
C-function counts degrees of freedom, at least near the 
critical points where it equals the central charge of the 
conformal (scale-invariant) fleld theory. We recall some 
elementary properties of the spectral form. In two di- 
mensions symmetries of space-time reduces the spectral 
function p to a single scalar function, which by causality 
(unitarity) is positive definite. This function measures 
the density of degrees of freedom in the given theory, 
with poles at single-particle states and cuts at the con- 
tinuum, which sum to unity. If the theory is changed 
so that a new set of low-energy states become relevant, 
then this theory has a different spectral function. We can 
imagine that these two theories belong to a single, one- 
parameter family of effective actions by a suitable choice 
of basis (lagrangian parameters or coupling constants), 
giving a family of spectral functions S R (a S R) that 
interpolates smoothly between the two theories. The cor- 
responding family of RG potentials $a is obtained by in- 
tegrating the spectral function. Friedan's dissection of 
this potential shows that it must be a positive real func- 
tion, and the C-theorem shows that the RG flow must 
be a gradient flow, i.e., "downhill". As long as the ef- 
fective theory is well behaved this should be true in any 
dimension since it avoids limit cycles and other patholo- 
gies of the flow. The existence of these potentials is a very 
useful simplification, since it lifts the geometric analysis 
of renormalization from vector to scalar fields, and RG 
potentials will therefore be the starting point for our gen- 
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eralized analysis presented below. 

To this reasonable list of properties we add that the RG 
potential should be "holomorphically factorized", which 
means that: 



$((7, a) = In |(y5p = In (y5((T) + \nip{(7) 



(3) 



It is really the partition function Z(cr, ct) = |CP = 
C(cr)C((T) that factorizes into a holomorphic and anti- 
holomorphic part, but the terminology is inherited by 
the "thermodynamic" potentials derived from InZ. The 
vacuum energy, for example, is F{a,a) oc lnZ(cr, (f) = 
ln(^((T) + lnC((T), and we expect the "vacuum entropy" 
<i>(cr, (f) to "factorize" in the same way. This is automat- 
ically the case at quantum critical points (because of the 
local conformal symmetry), and it also agrees with the 
flows observed in spin-polarized quantum Hall systems 
where a holomorphically factorized potential appears to 
account for all available scaling data. 

It now follows from eq. ([2| that the covariant (3- 
function should be the logarithmic derivative of a holo- 
morphic potential ip: 



Pa 
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(4) 



For a system with modular symmetry this implies that 
$a is a real-valued modular function that is finite in the 
finite part of H, and fi^ is a holomorphic weight 2 mod- 
ular form. A modular function is not, strictly speak- 
ing, a weightless modular form, since it must have a 
pole somewhere, but since we are studying the system 
at strong coupling it is physically reasonable to push this 
singularity to infinity, where the model decouples and we 
expect the cr- model to be irrelevant. The existence of 
these potentials is automatic in modular mathematics, 
since weight w — 2 forms actually do derive from holo- 
morphic potentials, showing again how remarkably well 
suited modular mathematics is for the task at hand. 

In short, we have argued that a physical RG potential 
is given by a holomorphic modular function ^p{cr), as in 
eq. ([3|. While the argument was conceived in the con- 
text of the spin-polarized QHE, as presented here there 
appears to be no obstruction, physical or mathematical, 
to considering more general situations. It is therefore 
our purpose to test this conjecture in the widest avail- 
able context. In other words, emboldened by the suc- 
cess in the spin-polarized case, we seek other situations 
where the effective theory may be harnessed by an emer- 
gent modular symmetry. The first step is obviously to 
consider other quantum Hall experiments, where addi- 
tional degrees of freedom are relevant. This includes sit- 
uations where spin is not fully polarized, multi-layered 
systems, graphene, etc., collectively referred to as "multi- 
component quantum Hall systems" in the following. In 
order to develop the vocabulary needed to discuss their 
symmetries it is advantageous to first summarize some 
basic facts about parameter space symmetries, modular 



transformations and RG flow geometry in the next sec- 
tion. Equipped with the appropriate mathematics we 
then proceed to a detailed symmetry analysis of the ob- 
served phenomenology of some multi-component systems 
in Sect. 3. 



II. SYMMETRIES AND RG POTENTIALS 

Mathematically, the full modular group r(l) = PSL(2, Z) 
is the simplest of the modular symmetries, but it does not 
agree with the scaling data in quantum Hall systems. 
It is, however, the starting point for the conventional 
mathematical discussion of these symmetries, reviewed 
here, which systematically develops the theory of mod- 
ular transformations from the group theory of r(l) and 
its subgroups. This may seem unnatural from a physics 
point of view, since r(l) is not a physical symmetry, but 
it gives a clear picture of how these symmetries are re- 
lated. This is crucial in our work, since we are looking for 
a framework that interpolates between the different effec- 
tive theories that are relevant for systems with different 
low energy degrees of freedom. 

We are fortunate that modular symmetries, i.e., sub- 
groups r C r(l), appear prominently in many branches 
of mathematics, and therefore are extremely well stud- 
ied. Since there are infinitely many modular subgroups, 
it is also very fortunate that scaling experiments severely 
constrain the physically relevant symmetries. As already 
mentioned T cannot be too large, but it also cannot be 
too small, since it then does not have natural candidates 
for the quantum critical points of the plateau-transitions. 
We can therefore confine attention to the most symmetric 
situations that yield phenomenologically viable RG flows, 
first considered systematically in the context of quantum 
phase transitions in Ref.[T31 This includes only "level 2" 
groups, which are the largest subgroups of r(l), as we 
now explain. 



A. Duality 

r(l) = (T, S) is generated by translations T{a) = cr + I 
and the simplest duality transformation S{a) = — f/cr, 
which satisfy two (and only two) algebraic constraints: 
5^ = 1 = [STY . They can be composed to give any inte- 
ger value of the matrix elements in 7 = (a, b; c, d) e r(l). 
This symmetry turns out to be too strong for physical 
applications, at least the ones considered here, so we con- 
sider the largest subgroups. They are obtained either by 
weakening the translation symmetry to T"(cr) = a + n 
for some integer n > 1, or by relaxing the duality sym- 
metry, or both. 

By a duality D we mean a transformation analogous 
to Kramers- Wannier duality, so that acting with D twice 
we get back to the starting point {D^ — 1). Conju- 
gating S with any group element X gives a transfor- 
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mation Dx = XSX ^ that is self-dual because S is 
(5- = S--!). The other constraint, {STS){TST) = 1, 
suggests that R — STS (the S'-conjugate of T), or equiv- 
alently its inverse W — R^^ = TST, will figure promi- 
nently together with S and T. On the upper half plane 
i?"(o') = (7/(1 — na). The main congruence subgroup at 
level two is r(2) = (T^,i?^), which means that the ma- 
trix representation only contains matrices that reduce to 
the identity mod 2 (i.e., a,d — 1; b,c — mod 2). 

There are four groups of modular symmetries "be- 
tween" r(l) and r(2), i.e., sub-groups of r(l) containing 
r(2) as a sub-group, which are called "level 2" symme- 
tries. The largest of these (it has index 2 in r(l), see 
for example Ref.[T3] [p. 5] for a definition) is Tp = r2 = 
{ST,TS). It is too large for our physical applications 
since it does not admit a physical potential. If it did the 
gradient of this potential would be a weight w = 2 mod- 
ular form, and no such form exists for this group. r(l), 
which contains this sub-group, is eliminated from further 
consideration by the same argument. 

If we want to keep all translations generated by T, then 
the group is unique: Tt = (T, Dr) = (T, i?^). This is es- 
sentially (up to outer automorphisms) the spin-polarized 
quantum Hall group. Keeping instead the original du- 
ality S requires that we double the translations to T^, 
giving the group Ts — (T^ , S) . Finally, the last of the 
groups at this level also has less translation symmetry, 
Tij = {Ry Dt) = {T'^jR). These three maximal symme- 
tries are viable candidates for physical symmetries be- 
cause they have physical potentials. We shall in fact 
see that these potentials belong to a r(2)-invariant one- 
parameter family that interpolates between these points 
of enhanced symmetry. 

If one of these symmetries emerges at low energy for a 
given physical situation, it is most easily recognized by 
identifying the parities of the attractors on the real line, 
i.e., the parities of the plateaux values ® — p/q. The 
modular group r(l) does not distinguish between the 
parities of the fractions p/q € Q, so all rational numbers 
are equivalent under this symmetry. r(2), on the other 
hand, respects the parities of both p and q, so it parti- 
tions the rationals into three equivalence classes. Each of 
the index 3 groups partition the rationals into two equiv- 
alence classes. With "o" representing odd integers and 
"e" representing even integers: 

rT = ro(2)-(T,i?2): {qee}U{qeo} 

rfl = r"(2) = (T2,i?) : {pee}U{peo} 

Ts = re(2) = (T2, S) : {pq e e} U {pq G 0} , 

where we have included synonymous group names fa- 
vored by mathematicians. 

A symmetry F identifies points in H, tesselating it with 
copies of the fundamental domain Fr, which can be cho- 
sen to be any subset of H where distinct points cannot be 
connected with transformations in F. The larger the sym- 
metry the smaller the fundamental region is, and if F has 
index /i in r(l) then Fr contains fj, copies of Fi = Fr(i'). 
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Figure 2; (Color online) Fundamental domains Fi and F2 
for the main congruence groups r(l) = PSL(2,Z) and F(2), 
showing conifold points of order two (® G H) and three (® G 
H), and the real cusps (0 G R). Not shown is the inequivalent 
cusp at infinity (ioo = 0). 



Fundamental domains Fi and F2 for the main congruence 
groups F(l) = PSL(2,Z) and F(2) are shown in Fig.[2] 
and fundamental domains for congruence subgroups of 
index /i = 3 are shown in Fig.jSj 

The boundary of a fundamental domain may contain 
conifold points where the symmetry has "folded" a local 
disc into a conical shape. When this occurs the fun- 
damental domain is a conifold, i.e., a generalized mani- 
fold with conical singularities that look locally like cones. 
This happens at points that are fixed by some group ele- 
ment X, and the deficit angle of the cone is (1 — 1/n) 27r, 
where n is the order of the group element, X" = 1. 

Fi has two such points. Since i = = (g) is fixed by 
the order 2 group element S (S*^ = 1), the deficit angle 
is TT. Similarly, since j = exp(27ri/3) = ® is fixed by the 
order 3 group element ST {{ST)'^ = 1)? the deficit angle 
is in total Att/3. 

We will always choose a fundamental domain that con- 
tains "the point at infinity". This point ioo, and possibly 
some of its real images, are called cusps (©), and by def- 
inition modular forms must not grow too fast near these 
points. IR stable (®) and UV unstable fixed points (q) 
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Figure 3: (Color online) Fundamental domain (a) F_r for Tr 
(b) Fs for Fs (c) Ft for Ft. Conifold points ® G H and real 
cusps are exhibited, the inequivalent cusp ioo = is not. 



of the RG flow are cusps. 

The pale critical points shown in Fig. [2] are excluded 
from the fundamental domains to avoid double counting: 
(a) j + 1 ~ j = ® since T e r(l), (b) 2 ~ = 
since G r(2). The dashed parts of the boundary 
are excluded for the same reason. Consider for example 
the solid semicircular boundary arc between and 1 in 
Fig.[2jb), given parametrically by z{9) = (1 + exp(i0))/2 
with 6 S (0,7r). Using the product of the generators 
of r(2) this maps to T'^R^{z) = (2 - 3z)/(l - 2z) = 
(3 — exp(— j0))/2, which parametrizes the dashed semi- 
circular arc between 1 and 2. Since z = 1 is the fixed 
point of T^B?, the dashed arc is a mirror image of the 
solid arc, both originating at this point. 

Regions enclosed by dashed and solid lines are images 



of the hyperbolic triangle A{i,j,oo). For aesthetic rea- 
sons we have displayed Fi = A{i + oo) U A(i -I- 
l,j + l,oo), rather than the more conventional choice 
A(i, j — 1, oo) U A(i, j, oo). F2 is six times as big (12 tri- 
angles) as Fi (2 triangles) because r(2) has index — 6 
in r(l). In the hyperbolic (Poincare) metric all trian- 
gles have the same area, but clearly not in the Euclidean 
metric used in these diagrams 

Similarly, the dashed boundaries and pale critical 
points shown in Fig. [3] are excluded from the fundamen- 
tal domain to avoid double counting: (a) 2 ~ = 
since e Fj^, (b) z + 2 ~ i = since e F5, 
(c) (3 + i)/2 ~ {1 + i)/2 = since T e Ft. F^ , ¥3 
and Ft are three times as big (6 triangles) as Fi because 
Tfj , Ts and F^ have index 3 in F(l). 

The subgroups considered here can only have conifold 
points of order 2. We shall see that /3-functions must 
vanish at these points, but the converse is not true, i.e., 
/3-functions can vanish at points that are not fixed by 
the group. Only in maximally symmetric cases (/x = 3) 
are RG and F fixed points the same. We use the icon 
to represent these critical points, whether they are fixed 
points of the group or not. 

Each symmetry F C F(l) determines a unique pair 
of phase diagrams, depending on whether the cuspoidal 
fixed point ioo ~ 1/0 = is attractive or repulsive. If F 
acts on the parameter space of a model, the space of con- 
ductivities, say, then this is the repulsive UV fixed point 
{ioo = 0) around which a perturbative expansion is usu- 
ally developed. But if F is acting on the inverse transport 
tensor (the resistivity), then it may be an attractive IR 
fixed point [ioo ~ 0). Given this one bit of information 
F determines the topology of the phase diagram. This 
is particularly easy to see in the maximally symmetric 
cases where F = Tx = R, S, T) . The phase diagram 
is dictated by which phases and which phase transitions 
are possible. If ioo is an attractive (repulsive) fixed point, 
then a transition 



/ ^p/q 



exists for: 



Ft ifi^ q and q' are even (odd) and 6 = 2(1), 

Fr iff p and p' are odd (even) and S — 1(2), 

Fs iflt pq and p'q' are even (odd) and (5 = 1(2), 

where S = det(/'; /) — det{p' ,p; q' , q) = p'q — pq' . These 
"diagnostic rules" follow immediately from the parity 
properties of the group, so the symmetry determines, and 
is determined by, the fixed point structure. 

Since the index 3 symmetry groups F(2) C Fx C F(l) 
[X = R, S, T) all have physical potentials they are the 
focus of our phenomenological analysis. These poten- 
tials belong to a F(2)-invariant one-parameter family, 
constructed next. 
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B. Analytic structure 



other singular points p are in the finite part of F, and 



We have so far concluded that we need to investigate 
the existence and uniqueness of mcromorphic potentials 
and holomorphic vector fields that are automorphic un- 
der some modular symmetry. The structure of a mcro- 
morphic function f{z) is completely determined by the 
points where it is not finite, i.e., its zeros and poles. By 
the order Up of a point Zp in H we mean the leading (low- 
est) power of the Laurent expansion at that point, i.e., 
f{z) oc (z — Zp)"p U Hp = then f{zp) is finite 

0, oo) and Zp is a regular point. If Hp is positive (and 
finite) then f{zp) vanishes and Zp IS 3t simple (rip = 1) 
or multiple (rip > 1) zero. If Up is negative (and finite) 
then f{zp) is singular and Zp is a simple (rip = —1) or 
multiple (rip < — 1) pole. If rip = ±oo the function is not 
mcromorphic. 

For functions that transform with weight w under a 
modular group the order of points on the boundary of 
the fundamental domain is not always obvious. Consider 
first the cusp ioo, which is the only cusp for r(l). We 
cannot Laurent expand f{z) in z— ioo, so by the order riao 
we mean in stead the leading power of an expansion in a 
suitable finite local variable. Because level N subgroups 
contain , modular forms at level N are periodic with 
period 27r7V, and can therefore be Fourier expanded in 
Qj^ = exp(2'Ki/N). The order rioo is the leading order in 
this qjy-expansion. 

For iV > 1 there are additional cusps on the real 
line that are images of the cusp at infinity. The order 
ric of a weight w function f{z) at a cusp Zc = 7c(«oo) 
is defined as the leading order of the g^y-expansion of 
{cz+d)~'^ /{"jciz)). For example, at level N = 2 the order 
no of f{z) at the origin — S'(ioo) is the leading power 
of z^™/(— l/z) expanded in powers of q2 — exp(i7rz). 

Cauchy's theorem can be applied to a mcromorphic 
function that transforms with weight w under a sub- 
group of index /i. By choosing a contour along the edge 
of F, taking care to circumnavigate the cusps and coni- 
fold points counter clockwise and recalling that conifold 
points contribute less than a full arc, it follows from 
Cauchy's theorem thatP^ 



l2 ' 



(5) 



where np is the contribution from the finite part of H. 

is the contribution from special points (that are fixed 
by elements of the group) on the boundary of the funda- 
mental domain, 



is the total order of these points. The contribution from 
Up will typically be severly constrained, or even elim- 
inated, by the physical requirements of RG potentials 
and flows. 

This sum rule is very useful since it allows us to classify 
the possible analytic behaviour of scalar and vector fields 
(potentials and /3- functions), i.e., the fixed points and 
singularities of effective field theories with F-symmetry. 

For the full modular group r(l), the contributions to 
the sum rule are 



2 



12' 



Subgroups at level 2 of index 3 have two inequivalent 
cusps, one at ioo and another at either or 1, and one 
conifold point of order two. The sum rule is therefore 



np = - (* = 0,1) 



For low values of w this sum rule is quite difficult to 
satisfy. For example, since the /3-function has weight 2 
and the cusp orders and Hp are integers, we must have a 
critical point at ig) (i.e., a simple zero contributing = 
1). Furthermore, since this function should derive from 
a physical potential we are led to consider the simplest 
possible analytic structure that saturates the sum rule. 



(no, no 



,) = (0,0,1,0). Remarkably, such entire 



and essentially unique /3-functions do exist, as we shall 
see below. 

The only subgroup at level 2 of index 6 is the main con- 
gruence group r(2), which has three inequivalent cusps 
at 0, 1 and zoo, but no conifold points. In this case the 
sum rule is 



no -I- ni 



w 
2 



which is a little easier to satisfy than the index 3 case. For 
example, there are now two possible entire /3-functions, 
which either has a simple zero at one of the cusps and no 
other zeros, or a simple zero in the finite part of F and 
no other zeros. 



2 3 



where n^ is the total order of all conifold points of order 
2, and n^ is the total order of all conifold points of order 
3 (® contributes iff the symmetry is r(l) or Tp). Any 



C. A potential family 

The modular groups we have been discussiii^ 
r(l), Tp, Tp, Ts, Tp, r(2), have canonical potentialaH 



9 





a 


TO 




^iTO/12 


(7 = 


no 


cr = 1 




o — ioo 




a = i 


Hi 




^3,3' 


Ub 


Up 


® G F 


J 







1 





OO 


-1 


oo 


-1 


oo 


-1 


1 








3 

















2 





1 





1 





1 





-1 





0, oo 


±3 












-1 





3 





oo 


-2 





1 





1 


2 





-1 











l + i 


Vs 


2 





3 








1 


oo 


-2 





1 


-1/4 





-1 











i 




1/2 





3 








1 





1 


oo 


-2 


-2 





1 











(1 + 0/2 


A 


1 





6 





1 





oo 


-1 





1 


1/2 


















fa 


a 





6 










oo 









* 




* 








(a) 


et 


— oo 


2 


6 


1 


* 





* 








1 


* 





* 





1 







et 





2 


6 


1 


* 








1 


1 





* 





* 





1 







et 


1 


2 


6 


1 





1 


* 





1 





* 





* 





1 







Pa 


a 


2 


6 


1 


(1 - a)* 





a * 





1 





* 





* 








1 


0-® (a) 


Pr 


-1 


2 


3 


1/2 


* 





* 





1 





* 





* 





1/2 





l + i 


Ps 


2 


2 


3 


1/2 


* 





* 





1 








1 


* 





1/2 





i 


Pt 


1/2 


2 


3 


1/2 


* 





* 





1 





* 





* 





1/2 





(l + j)/2 



Table I: Analytic structure of some modular functions and forms at level 2. A star (*) represents a "finite" but uninteresting 
value of the appropriate uniformizing variable. Pa = ii^ifit ~ "-^t) is "physical" iff o 7^ 0, l,oo, in which case it has a critical 
point that is on 9F if a is real. The last column contains the position of the saddle point of the potential, or equivalently, the 
location of the simple zero of the /3-function derived from this potential. J has a double pole at ioo and triple zeros at both 



the conifold points j and j' = 
thus saturating the sum rule. 



thus saturating the sum rule. Similarly, ipp has a triple zero at j and a triple pole at j' , 



(also called invariant or automorphic functions) 



^1 = (a2-a + i)3a-2(a-i)-2 

VP = (A + jf(A-/)-3 

= A(A-l)-2 
^5 = A(A-l) 

^2 = A = (02/03)' , 



(6) 



by their Fourier expansions: 

^i{u,q) = 291/4 ^(-l)"g"("+i)sin((2n + l)w) 

n=0 

00 

^2{u,q) = 2gi/'*^g"("+i)cos((2n + l)u) 

n=0 
00 

i?3(u,9) = 1 + 2 ^ g"' cos(2nu) 

n=l 
00 

^4{u,q) = l + 2^(-l)"g"'cos(2nu) . 



where 



The analytic structure of these potentials, together with 
some other functions that will make an appearance be- 
low, are listed in Tab.jlj 

The r(l)-invariant function ip-^ is Klein's famous J- 
invariant, J — (4/27) ipi . It has been expressed in terms 
of the r(2)-invariant function A, which is the ratio of the 
Jacobi ^-functions 02(0") = ^2(0, q) and Osier) = 1^3(0, g). 
The elliptic i?-functions -di are most conveniently given 



Since q = exp(7ricr) = exp(7ria;) exp(— Try), these sums 
converge very rapidly for reasonable values oiy — Scr and 
are therefore well suited for numerical work, especially 
graphics. 

J oc ip-^ and ipp are too singular to give rise to a mod- 
ular form of weight w — 2. There are in fact no weight 2 
forms at all with the symmetries of these potentials, so 
both r(l) and its largest (index fi ^ 2) subgroup Tp are 
too large for physical applications and we need consider 
them no further. 

The index three (/i = 3) subgroups are just barely 
small enough to admit weight 2 forms, but only one each, 
with a single critical points (simple zero) that is forced 
to coincide with a fixed point of the group. An RG flow 
diagram with symmetry Tn, Ts or Tt is therefore ex- 
tremely rigid, predicting the exact location of all critical 
points. The index six (p = 6) potentials are more flexi- 
ble, since there are two independent ones. They give rise 
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to a 2-dimensional space of weight 2 forms, which we can 
choose to span with two of the Jacobi 6'-functions, say 6*1 
and Thus it follows from modular symmetry alone 
that 9f must be a linear combination of these, a fact dis- 
covered by Jacobi and recorded in his famous identity 

The purpose of all this mathematical machinery should 
now becoming clear. Of the six possible potentials at 
level two, and Lpp are unable to satisfy the physical 
constraints placed on their analytic structure. The sim- 
plicity and similarity of the remaining four potentials is 
striking, and suggests an obvious way to embed them 
into a single, one-parameter family of physical poten- 
tials. Consider first the general r(2)-invariant RG po- 
tential — ip -\- ij: obtained by a linear superposition of 
the maximally symmetric potentials in eq. (|6| : 

V'(o') = In ip^ -\- Cglnipg + Crphi ifrp , Cj^ g J, C 

This family interpolates smoothly between the maxi- 
mally symmetric cases, but there are too many param- 
eters for it too be useful as it standJ^^. We can ignore 
an arbitrary normalization factor that only affects flow 
rates, not their shapes, leaving two (complex) parame- 
ters. However, since (/?^, ipg and ipj, are simple rational 
functions in the fundamental r(2) modular function A, 
the level two symmetric family of quantum Hall RG po- 
tentials reduces to (we choose a convenient parametriza- 
tion for a): 



$a oc lnA(A - 1)"-^ -hc.c. 



D. Modular RG flows 



(7) 



From eq. Q it immediately follows that the /3-function 
1 - aA A' 



/3a,, 



-9„<I>a(cr) OC 



A - 1 A 



inO^il - aA) (8) 



by construction satisfies all physical requirements, thus 
proving the existence of a one-parameter family of phys- 
ically viable /3-functions that interpolates between the 
enhanced symmetries T^^^.s- The second equality in 
eq. (|8]), which follows from the Jacobi identity and the 
fact that A' /A = inOf, is particularly useful for locat- 
ing the critical point. Since 0- functions are finite on F 
(comp. Tab.|l|, (3 can only vanish if A(cr0) = 1/a, which 
can be inverted in terms of elliptic integrals of the first 
kind. Consequently, the /9-function has a simple zero at 



= A-i(l/a) 



iK'{l/a) 
K{l/a) 



mod r(2) , (9) 



as it should, and nowhere else (in the fundamental do- 
main), as it should. 

The r(2) symmetry of $a is enhanced when a takes 
certain values (Fig.|4]): to Tji when a = aj^ = —1, to Tx 
when a = ot = 1/2 and to when a — as — 2. The 
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Figure 4: The flow of the vector field Pa is r(2) symmetric 
for any value of the parameter a G C. When a = —1, 1/2, 2 
the symmetry is enhanced to an index 3 group, and when a = 
0,1,00 the /3-function degenerates to a ^-function, which is 
finite on the interior of the fundamental domain and therefore 
has no quantum critical point. 



corresponding /3-functions reduce to particularly simple 
Eisenstein functions^: 



oc 



I3s oc 



(3t oc 



A + 


1 


A - 


1 


2A- 


- 1 


A- 


1 


A- 


2 


A- 


1 



oc 1 + 



nq 



oc 1 



oc 1 



f. (2n + l)g^"+i 



nq 



2n 



^2T^ 



n=0 

has a Fourier expansion in q^ — exp(27rif7) because 
T £ Tt, while f3fi and f3s are Fourier expanded in q 
because the smallest translation in F/j and Tg is T^. 

When a = 0, 1, oo the /3-function degenerates to one of 
the Jacobi 6'-functions, 

Pa=o OC = ^n9t 



Pa=l OC 

/3a^ 



A 



— > ITT I 



Since these are finite on the interior of the fundamen- 
tal domain the would-be quantum critical point has 
been pushed off the physical (finite) part of parameter 
space. The uniqueness of the interpolating potential 
in eq. ([?]) follows from the physical requirements we place 
on the analytic structure of the /3-function. We want a 
simple zero at the delocalization critical point a(g, at a 
finite value of 5(t, and the mildest possible singularities 
allowed by modular symmetry. Rankin's theoremJSI ^gU 
us that the weight of (3^ can be carried by the logarithmic 
derivative 9^ In A = A'/A = inOf of the level two invariant 
A, and that the remaining weightless factor is a rational 
function of that invariant. By inspection of /3fl, f3g and 
/3t it is clear that a fractional linear form is needed, but 
no more in order to avoid non-simple critical points. In 
other words, /^^ oc F{\)dcr In A oc 6fF{X) with 



F{X) 



AX + B 
CX + D 



Since the weight 2 form 6f is finite (no zeros or poles) 
on H, any non-trivial analytic structure is carried by F. 



11 



The modular invariant A has a simple pole in q on the 
fundamental domain F2. This cannot be avoided since 
the only holomorphic invariant is a constant. We now 
use the sum rule eq. ([s]) for modular forms that transform 
with weight -u; = or w = 2 under r(2) = 6): 



no + ni + rioo + 



B/A 



The possible zeros and poles of F{X) are at A 
and A = —D/C. If the "particle-hole" symmetry /3(A) — 
/3(A) holds, then A, B, C, D are restricted to be real, and 
zeros and poles can only appear along the semi-circles 
where A is reaP^. The requirement of having no poles 
except at the boundary IHI^ sets —D/C = 0, 1 or 00. 
Similarly, the requirement of having a zero on H, sets 
-B/A 0,1,00. If D/C = -1 there is a pole in F 
at a = 0, so no{F) — —1. However, this does not af- 
fect the analytic structure since Of has a zero at cr = 0, 
so 71-0(^4) = 1- The pole from F is therefore cancelled, 
no(/3o-) = 0, and this type of /3-function can have the 
required critical point. 

On the other hand, if D/C — there is a pole in F 
at cr = ioo, giving nao{F) = — 1, but since nooiOf) — 
this contribution is not cancelled, n-ooiPa) = —1, and 
there must be additional zeros to saturate the sum rule 
for the /3-function. Since we only want one simple zero, 
this is not a viable form of the /3-function. Similarly, if 
D/C = —00 then ni{F) = —1, and since ni{9f) = 0, 
this also requires additional zeroes and is therefore not a 
viable form for a physical ^-function. 

Since the metallic conductivity point (cr -> ioo) should 
be a repulsive fixed point, the /3-function has the form: 
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A' 



A(l-A) 



{AX/B + 1) = -in9j{AX/B + 1) , 



up to an undetermined positive multiplicative constant, 
which allows us to rescale A/B = —a, giving eq. ([8|. This 
fixes /3a to be the gradient of the ansatz potential exhib- 
ited in eq. Q. 

When a =/= {aR,aT,as} = { — 1,1/2,2} the flow is 
only r(2) symmetric, and the critical point with finite 
imaginary part (Qa > 0) is not a fixed point of the 
group. However, because A is real along the boundary 
of ¥fSi (p. 228), as long as a S M the critical point is 
still on the semi-circular boundary of F2. The semi- 
circle law observed in the spin-polarized case, for which 
a = ay = 1/2 S R, is therefore an automatic consequence 
of modular symmetry, as pointed out in Ref. 15, 

Fig. [5] shows a sequence of ten distinct RG fiows, ob- 
tained for real values of the parameter in the range 
a G [—10,-1-10], with the choice ioo = Q (if ioo — ® 
the direction of all arrows is reversed). The location of 
the quantum critical point (X) (and some of its modular re- 
fiections) is computed using eq. ^ , and as expected they 
move along boundaries of fundamental domains for r(2). 
Thick (blue) lines are phase boundaries or separatrices. 




Figure 5: (Color online) r(2)-symmetric flows generated by 
for ten values of a £ [-10, 10] (e = 0.05). At a = -1, 1/2 
and 2 the symmetry is enhanced to Tr, Ft and Fg. Thick 
blue lines are phase boundaries or separatrices. 
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For increasing values of a < the quantum critical 
point ® first moves straight down the vertical line SRcr = 
1 (a = — 10, — 1, — £ = —0.05), until it fuses with the 
repulsive point = lasa— >-0~, denoted by ® = 1 in 
the a = diagram. A further increase in a appears to 
split ® into three parts {a — e): a mirror pair of critical 
points that move apart along semi-circles (as described 
in detail in Sect. 2.1: 0+ ^ 0, T'^R^{®) 2"), leaving 
behind a new attractive point © = 1. This attractor is 
attached to a new phase injected during the split, which 
squeezes out the old phases (a = 1 — e) as a increases 
and eventually engulfs much of the diagram (a = 1). 
At a = 1/2 this new phase provides the "real estate" 
needed to enhance the translation symmetry from to 
T (and the total symmetry from r(2) to Tx)- When the 
parameter reaches the value a = 1 the exact opposite of 
a split occurs: mirror pairs of critical points coalesce at 
® = 0, and then move straight up (a — 2, 10), leaving 
behind a repulsive fixed point = 0. 

Notice how this resolves the enigma of having a discon- 
tinuous "jump" in symmetry ata = — 1, 1/2 and 2, which 
at first sight seems mysterious and unphysical, especially 
since the model acquires an infinite number of new non- 
commuting symmetries. However, our discussion of the 
RG potential $a, which generated the flows shown in 
Fig. [5] makes it transparent that this "spontaneous sym- 
metry generation" is no more mysterious, discontinuous 
or "spontaneous" than ordinary "spontaneous symme- 
try breaking" . The potential changes continuously and 
smoothly as a function of a, interpolating between differ- 
ent (infinite and non-abelian) symmetries by switching 
parts of the potential on or off, thus injecting, expanding 
or shrinking phases as needed^'^. 

Some of the diagrams in Fig. [5] are similar to those 
appearing in Ref.l45[ where a different interpolation be- 
tween enhanced symmetries based on "scaling functions" 
is investigated. We defer a more detailed comparison to 
Sect.HVl 



III. COMPARISON WITH EXPERIMENT 

We wish to confront the proposed RG potential <i>a 
(eq. ([7|) with available scaling data. The simplest cases, 
when the spin is fully polarized or fully resolved, have 
previously been identified as having the e nhanced sym- 
metries Ft (at = 1/2) and F;^ [a^ = -l^pHMUl. The 
RG flows generated by $1/2 and $_i were found to be in 
excellent agreement with availa ble sca ling data for both 
integer and fractional transitioni^^'^ . 

Here we extend this analysis to more general situations 
that allow us to test if the proposed potential does inter- 
polate between the points of enhanced symmetry consid- 
ered previously. Since we shall be discussing both flows 
in the complex conductivity parameter a and resistivity 
parameter p = S{a), when confusion can arise the critical 
point ® and the parameter a will be indexed accordingly. 
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Figure 6: (Color online) Temperature driven flow datsP^ for 
(a) the ®CT 1 transition, and (b) the 1 (8)0- 2 transition, in an 
AlGaAs/GaAs heterostructure with density n — 2x 10^^/cm^ 
and mobility fi — 7.6 x lO^cm^ /Vs. The insets show the loca- 
tion of the quantum critical points in the conductivity phase 
diagram, compared with the critical points of the interpolat- 
ing RG potential $0.89. 



We are not at this time able to derive the value of 
a from micro-physical considerations. The question we 
address here is whether it is possible to parametrize all, or 
almost all, scaling data with a single real number a that 
fixes the modular RG potential $0 , and therefore the 
location of all quantum critical points and the geometry 
of all RG flow lines. This "phenomenological" approach, 
explained at length in the introduction, is extremely rigid 
and therefore eminently falsifiable. 



A. Spin split Landau level 

Huang et al. have studied the conductivity transition 
between the first and second Hall level, which we denote 
by 1 2, and also the (So- 1 transition, in a semi- 
conductor sample where the spin degeneracy of Landau 
levels is small but resolved. The gap between the two 
spin flavors is determined by the ratio of the electron ef- 
fective mass and effective g* -factor in the sample. The 
sample is a AlGaAs/GaAs heterostructure with density 
n = 2 X 10^^ /cm^ and mobihty = 7.6 x lO'^cm'^ /Vs. 
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For the (X)p 1 transition they reported a a critical 
magnetic field value ~ 14.7 T, a critical value of the 
diagonal resistivity p£i(0 (g) 1) « 37fcr2 « 1.43 [h/e^], and 
that ph appears to be independent of B for the lowest 
temperatures (see Fig.[6]ja)). 

However, the crossing point appears to be moving 
down with decreasing temperature, and we estimate that 
-Bx (0 ® 1) w 14.6 ± 0.1 r, which gives the critical value 
pg, sa 1.007 + 1.439z [h/e^] for the Q®pl transition. This 
is very close to the value ~ 1 + 1.439i [h/e^] obtained 
from $a with ap ss —5.25. The corresponding value for 
the conductivity flow parameter is k, 0.84, and the 
phase diagram is shown as an inset in Fig.[6f^a). 

For the 1 (8)p 2 transition they find that the critical 
magnetic field value is By^ « 5.55 T. This gives a critical 
point on the semi-circle shown in the inset in Fig.[6][b), 
which is almost identical to the critical point of $o.84 ^-Iso 
shown in the same diagram. 

These results appear to show that there is a r(2) sym- 
metry emerging in this sample, with an RG potential 
$0.84 whose critical points are close to the observed val- 
ues for the 0(8)crl(8)o-2 transitions. Using the same sample 
as in the experiment described abov^i^l^ Huang et al.l^ 
have studied the temperature scaling in more detail. The 
experiment probed temperature driven RG flows, with 
temperatures sampled at T = 0.94, 0.68, 0.49, 0.31 [if], 
and their data are reproduced in Figs. [t] and [s] (adapted 
from Ref.JOj. 

The semi-circular geometry of the flows strongly sug- 
gests that the flow is following the separatrix connecting 
the plateaux at = 1 and = 2, if we allow for a 
systematic error that shifts ajj up by 0.17e^//i, in which 
case a real value of the parameter a will suffice. This 
value is most easily determined by finding the location 
of the quantum critical point cr^, and comparing this 
with the best fit to <i>a- If were known exactly we 
could simply compute the value of a using eq. ([9]), i.e., 
a = l/\{a^). However, with the substantial uncertain- 
ties in currently available data it is often better to make 
a global comparison of all the scaling data, not just those 
near the critical point, in order to find a reasonable value 
for a giving a good overall fit. 

For the 1 (g)cr 2 transition we find that a = 0.93 ± 0.01 
gives a good fit to the data, as shown in Fig. [7] where the 
fiow generated by <i>o.93 has been superimposed on the 
data from Ref.[2Ql Notice that there is a cluster of data 
that is almost stationary, and therefore should be very 
close to an RG fixed point. This is indeed the case, as 
the global data fit gives a quantum critical point (g) that 
almost totally echpses this data series. 

For the (8)o- 1 transition we find that a = 0.89 ± 0.04 
gives a reasonable fit, as shown in Fig. [8] where the fiow 
generated by $o.89 has been superimposed on the data 
from Ref. 20. Notice that there is also in this case a 
cluster of data that is more or less stationary, albeit not 
as convincingly as in the previous diagram. These data 
points are indeed quite close to the candidate critical 
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Figure 7: (Color online) Temperature driven flow data for 
the I® 2 conductivity transition (adapted from Ref. l20|l . com- 
pared with the theoretical flow generated by the interpolating 
RG potential $0.93 • Thick (blue) lines are theoretical phase 
boundaries. Dashed blue lines are separatrices. 



point (g) identified from the rest of the flow data, as ex- 
pected if the global fit is reasonable. 

In both diagrams the thin solid (red) fiow lines have 
been chosen to illustrate that the data are consistent with 
the potential. In Fig. [7] the fit is essentially perfect, even 
with the small error bars we have estimated (since none 
were published). Comparing these results with Fig.|6j 
there appears to be a r(2) symmetry emerging in this 
sample, with an RG potential consistent across sev- 
eral quantum Hall transitions. The small dispersion in a- 
values is not surprising given the precision of these exper- 
iments. Since we estimate that Bx (1 (g 2) « 5.54 ± .05 T, 
the uncertainty in Ua is ±0.1 in both transitions. 

In addition to the the fiows presented here, the tran- 
sition 2 (8)p 4 was also investigated in Ref.UHl They 
found a critical point at = (7.75, 3.54) ± 0.05 kil w 
(0.300,0.137) ± 0.002 [h/e^], which corresponds to Op = 
— 1.76. The authors remark that the plateau in pu is 
not well developed, and note that if we assume that the 
true value is po = 1/10 [h/e^], then pn — 3/10 [h/e^] (so 
p^ = {3 + i)/10 [/i/e^]) corresponds exactly to the pa- 
rameter value Up — —1, i.e., the enhanced symmetry 
for the usual spin unresolved quantum Hall effect, which 
is often observed at higher Landau levels. 
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Figure 8: (Color online) Temperature driven flow data for 
the 0® 1 conductivity transition (adapted from Ref. l20|l . com- 
pared with the theoretical flow generated by the interpolating 
RG potential $o.89- Thick (blue) lines are theoretical phase 
boundaries. Dashed (blue) lines are separatrices. 



B. Zeeman splitting 

Hang et al.'^ have studied the effect of Zeeman splitting 
in IQHE transitions in three different samples (called A, 
B and C) with various spin degeneracies in the Landau 
spectrum. Sample A is the same sample as in Refs. 1191201 
analyzed in the previous section. Sample _B is a GaAs 
system with InAs dots, but neither the mobility nor the 
density have been reported. Sample C is a SiGe hole 
system with mobility /i = 1.69 x Ky^cm? /V s and density 
n « 2.2 X 10^^/cm^. Their data are reproduced in Fig.^ 
While the data are less than ideal, they do appear to 
provide some evidence for a generalized semi-circle law, 
which asserts that all quantum critical points should ap- 
pear on the boundary of F. This is an automatic property 
of $0 flows if a is real. 

Only the plateau at « 2 is properly developed, 
while the second visible plateau at w 1 does not have 
a fully developed zero of po, indicating that the scaling 
regime has not been reached. Accordingly, the transi- 
tions deviate somewhat from the expected semi-circles, 
although they do appear to indicate the common trend, 
with two transitions above and two under the semi-circle. 
The critical points are on the semi-circles to a reasonable 
accuracy, within a couple of the rather large standard 
deviations in these experiments. 

Their conductivity data allows us to compare the lo- 
cation of the critical points in the transitions (8)p 1 <8)p 2 
for sample A, giving (8)p 1 w (0.31,0.44) ± 0.01 and 
1 (8)p 2 sa (0.75, 0.47) ± 0.01. These correspond fairly well 
to the value Oo- = 0.895 ± 0.025, in agreement with the 
estimates ~ 0.93 ± 0.01 and Ua = 0.84 ± 0.1 obtained 
for this sample in the previous section. 

The critical points for the other samples are also ap- 
proximately on the semi-circle. For sample C our esti- 
mates for the value of a are different for the 1 (g) 2 and 
2(8)3 transitions, but this is to be expected since the spin 



configuration drastically lowers the the splitting between 
higher Landau levels (see Ref.[21] Fig. 1). The 2(^3 tran- 
sition therefore only enjoys a r(2) symmetry, whereas 
the the well resolved transition 1 (g) 2 obeys the full Ft 
symmetry. 

For sample B similar remarks apply, but in this case 
the (E" 1 transition reflects the spin- split symmetry F^, 
while the 1(8)2 transition reflects the spin-unresolved sym- 
metry Fa (compare Ref.l^Tl Fig. 1). As the best fit to the 
data we find ~ 0.52±0.4 for the 0(^)1 transition. From 
the 1 2 data we extract w -0.94 w -0.94 + 0.46*. 
This means that the flow is very close to the a„ = —1 
flow near the semi-circle, and that the complex part of 
flo- only displaces the critical point and the flow slightly 
from the centre of the semi-circle. 

We interpret these experiments as providing support- 
ing evidence for the generalized semi-circle law. For 
sample A sufficient data are available for us to extract 
Oct = 0.895 ± 0.025 from the location of the critical 
points, and this is consistent with Ua ~ 0.93 ± 0.01 and 
Oct = 0.84 ±0.1, obtained in the previous section from 
the data in Refs. 1191201 
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Figure 9: (Color online) Critical points for the transitions 
studied in the samples A, B and C discussed in the text, are 
in rough agreement with the generalized semi-circle law (data 
adapted from Ref.[2T]Fig- 3). In order to make a global com- 
parison some of the data have been rescaled and translated 
to the proximity of the exhibited semi-circle. The margin of 
error is indicated by black crosses. 



C. Overlapping spin sub-bands 

Murzin et aX!^ have studied a quantum Hall system in 
which the spin splitting is small enough for spin split Lan- 
dau levels to overlap substantially. Their three samples 
(labelled S34, S40 and S50) are heavily Si-doped GaAs 
layers with mobilities /i ~ (2.0 — 2.4) x 10^ err? jV s and 
densities n = 3.9, 4.6, 5.0 x V^^^ jcrv? . They estimate the 
Zeeman splitting to be of the order of 4 K, whereas the 
level broadening is much larger, of the order of 100 K. 

Data obtained above a few hundred milli-Kelvin are 
not expected to be in the scaling domain, so we disregard 
the high temperature data (which start at T w 12 K). 
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Figure 10: (Color online) Low temperature flow data for the 
®CT 1 ®cT 2 transitions in sample S34 (adapted from Ref.l22l 
Fig. 2), compared with the theoretical flow generated by the 
interpolating RG potential <l?o.i8- Thick (blue) lines are theo- 
retical phase boundaries. Dashed (blue) lines are separatrices. 



The situation is quite similar to the original experi- 
ment where scaling was first observed in the QHE^"^ , and 
where this procedure was justified, but it has only re- 
cently been pointed out that scaling data exhibit a hyper- 
sensitivity to temperature that still does not receive suf- 
ficient attentiorP. 

The low-temperature data (down to T = 0.1 if) are 
reproduced in Fig. 10 (adapted from Ref. l^ . and com- 
pared to the flow generated by with a = 0.18 ± 0.01. 
Since the experimental errors (which are not reported) 
must be at least as large as the plot markers used in 
Ref."25', Fig. 2, there is no discrepancy between the data 
and the theoretical flow. 

The theoretical critical points ®cr 1 ~ 0.66 + 0.47i 
and 1 ®„ 2 ss 1.34 -t- 0.47* are reasonably close to the cor- 
responding experimental values uq^i = (0.5 ± 0.05, 0.5 ± 
0.02) [e^/h] and di^a = (1-5 ± 0.05,0.5 ± 0.02) [e^/h], 
based on a linear extrapolation of the (Tc-maxima to zero 
temperatur^22l_ 

In Ref.im it was suggested that their data could be 
interpreted as a hybrid flow composed of two spin com- 
ponents scaling independently. However, the good agree- 
ment we have found with the r(2) symmetric potential 
suggests that the observed scaling simply follows from 
reduced spin splitting that leads to near degenerate Lan- 
dau levels. The most obvious consequence of this re- 
interpretation is that the location of the quantum criti- 
cal points is shifted away from the more symmetric spin- 
polarized positions, as shown in Fig.[TOj This should be 
straight forward to verify by more accurate experiments. 



D. Layered heterostructure 



Huang et al.!^ have studied the Q ®p 1 ®p2 transitions 
in a layered AlGaAs/InGaAs/GaAs heterostructure with 
mobility /i = 1.0 x lO^crn^ /Vs and carrier density n = 



1.3 X 10^^ /cm?. Again the high-temperature data (above 
1 K) are not relevant for our study of the quantum phase 
transitions. 

For the 1 (gip 2 transition they observed a stationary 
flow that they interpreted as being controlled by a criti- 
cal point « 1. 9-1-0. 39i [e^/h]. In comparing this result 
with the RG potential $a, the best value we have found 
is a„ w 0.992, which gives w 1.85-f 0.35i. However, in 
this case the experimental flow does not fit the theoret- 
ical fiow. To avoid having experimental and theoretical 
fiow-lines crossing, the data might be conjectured to fall 
on a semi-circular separatrix, but then the "semi-circle" 
formed by the experimental data points would be skewed 
and distorted, so this seems implausible. 

For the (8)p 1 transition they report « 20 fcfi « 
0.77 [/i/e^], and (^p 1 « 1 0.77i[h/e'^]. This gives 
0(8)cr 1 ~ 0.62-f 0.48«, which corrresponds to a„ ~ —0.32, 
in contradiction with the value a sa 0.992 we obtained 
above for the 1 (g) 2 transition. The critical points of the 
two transitions appear at values of an that are too high 
for them to be related by r(2) symmetry, i.e, they are 
not symmetric with respect to the an — 1 line. 

The authors of Ref.[211 comment that these rather un- 
expected results could be due to a variety of reasons: 
inhomogeneity in the sample, mixing between Landau 
bands, small spin-splitting, or temperatures too high for 
scaling to be observed. In any case, this system does not 
appear to belong to the family of RG flows generated by 
^a, and therefore highlights the need for a better un- 
derstanding of necessary and sufficient conditions for the 
emergence of modular symmetries. 



E. Remarks on graphene 

A single layer (b — 1) of graphene has a quasi-particle 
spectrum with a "relativistic" dispersion relation, i.e., 
they are essentially Dirac fermions. Bilayer (6 — 2) 
graphene has the spectum of massive Dirac fermions. 
The free particle Landau problem gives integer quantum 
Hall (IQH) plateaux at: 



N eZ, 



(10) 



in natural units (e^/h). The factor 4 arises from two-fold 
spin and valley degeneracies in graphen^^D, 

As originally noted in Ref.[2Bl if a is replaced by <j/2 
in order to take account of either the spin- or valley 
degeneracy, then the fixed point structure of the QHE 
in single-layer graphene appears to respects the modu- 
lar symmetry Tg. The system can be argued to be a 
spin-polarized system, but because of a Berry phase it 
pick ed up by Dirac fermions moving around the Dirac 
con^^, the group Tt should be conjugated by a 7r-flux 
attachment STS, which gives the bosonic duality group 
p^2 13 27\^ This reduces the translation symmetry from T 
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to T^, and this leads to the correct periodicity of 4 in 



eq. (lOl 



In bilayer graphene the Berry phase is 27r, which gives 
the symmetry Tji after the conductance is divided by 2. 
This is consistent with eq. ( [To| , except that the cr© = 
state is missing. Instead there is a metaUic state with 
a minimum conductivity. The zeroth Landau level state 
can be created by doping. 

Recent advances in graphene sample fabrication have 
revealed fractional quantum Hall (FQH) states^^, as well 
as anomalous IQH states that do not fit into the non- 
interacting sequence given by eq. ( [lO| . 

The anomalous states have been related to symme- 
try breaking between the spin- and valley degeneracies 
of Dirac fermions at high magnetic fielda^S, and by now 
almost all IQH states have been observed in graphene de- 
vices at high magnetic fields-^ In particular, all the 
broken symmetry states of the iV = Landau level, giv- 
ing plateaux at cr© = 0, 1, 2, 3 and 4, have now b een ob- 
served for both single- and bilayer graphen d^^ l ^'^ l There 
is also evidence for all plateaux a© = 1, . . . , 16 for bi- 
layer graphend^H^ which suggests that the strong field has 
completely "polarized" the system, lifting the multicom- 
ponent degeneracy and returning the symmetry to the 
conventional one-component group Tt- Recently, Ref. l33l 
reported fractions that are predominantly multiples of 
1/3, although some states are missing from the = 
Landau level. 

Unfortunately scaling data for the plateaux transitions 
in graphene are not yet available. There is some evidence 
that the plateau-plateau transitions are in the same uni- 
versality class as in the usual spin-polarized QHE, re- 
ported in Refs. 1341351 but they also find evidence that 
the insulator-plateau transition is in a different univer- 
sality class. 

Nevertheless, the fixed points of the fiows appear 
to be consistent with the symmetry transmutation^'^ 

Tg — > driven by an increasing magnetic field, 

leading eventually to a single flavor of spin-polarized 
fermions without the 7r-Berry phase due to a (sponta- 
neous) mass gap. Clearly this kind of a scenario (or some- 
thing more complicated) can be achieved with the pro- 
posed ansatz, where a is a non-universal sample specific 
parameter that depends on impurities, sample chemistry, 
sample architecture, strength of the external magnetic 
field, Zeeman splitting, valley separation, etc. 

Finally we remark that if low energy models of 
graphene do enjoy modular symmetries, then they can 
in principle be used as a diagnostic tool to characterize 
the material. Abanin and Levito'v'^ll have recently inves- 
tigated another way to use transport measurements to 
characterize materials exhibiting the quantum Hall ef- 
fect, including the number of layers in a graphene sam- 
ple. Since the total macroscopic conductance measured 
in the two-terminal experiments considered by these au- 
thors depends on sample geometry (e.g., the aspect ra- 
tio) , it is conceivable that this will conceal the properties 



of the material (e.g., the electron density) and therefore 
render the method unsuitable for sample characteriza- 
tion. Using a semi-circle relation between the conduc- 
tivities they show that this is not the case, so that their 
method can be used to determine the number of layers 
in graphene, for example. 



IV. COMPARISON WITH OTHER WORK 

There have been various attempts to analyze the phase 
structure of the QHE, starting with a proposal by 
Khmel'nitzkiP^that identified a translation symmetry of 
RG fiows in the IQHE. This was motivated by a sigma- 
model of localizati on, w hich was extensively elaborated 
by Pruisken et al.'^SESl^ but the target space geometry 
does not appear to be rich enough to include the FQHE. 

Ref. 2] proposed that a modular symmetry would be 
capable of describing both the integer and the fractional 
Hall effects by including dualities in addition to transla- 
tions, as described in the introduction. The three max- 
imal index 3 subgroups Fx {X = R, S, T) of the full 
modular group PSL(2,Z) were immediately identified as 
the largest symmetries of physical interest. The group 
AutFT, acting by holomorphic fractional linear trans- 
formations on the complexified conductivity, was shown 
to give the correct phenomenology (fixed point structure 
and super-universality) of the spin polarized QHE. Since 
p = S{a) this is equivalent to having an emergent F^; 
symmetry (which is S'-conjugate to Tt) on the space of 
resistivities, and F5 was proposed as the relevant symme- 
try for analogous transport problems with bosonic quasi- 
particleP'. 

The symmetry F(2) and interpolations between the en- 
hanced symmetries Tx (X — R, S, T) were introduced in 
Ref.Hni where the idea of an interpolating potential elab- 
orated here was also proposed. 

At r oughl y the same time superficially similar 
dualities^SlEll acting on the filling fraction v were consid- 
ered, and the resulting transformations are known col- 
lectively as "the law of corresponding states" . Since v is 
essentially the plateau value u© = a^j £ M, this approach 
is oblivious to the complex structure that gives modular 
symmetries most of their predictive power. These duali- 
ties appear to disagree with experiment, unlike the com- 
plexified duality identified in Ref.|2 | that is in excellent 
agreement with available datcPE^EH, 

Following the ideas in Ref.'W, Burgess et al.'^S! re- 
lated complexified modular transformations to dualities 
in the low energy (linear response) transport coefficients 
of an effective Chern-Simons Landau-Ginzburg field the- 
ory. These dualities follow from periodicity of the statis- 
tics angle of the charge carriers in two dimensions, to- 
gether with particle vortex duality, assuming that only 
particles or only vortices take part in low energy charge 
transport. This gives symmetry for fermionic charge 
carriers, and F5 symmetry for bosonic carriers. Note that 
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dualities relate systems with different (anyonic) statistics 
of charge carriers, but never bosonic charge carriers with 
fermionic ones. These dualities were further extended to 
the non-linear regime of charge transport in Ref.B51 

Burgess et al.^ have also considered generalized mod- 
ular symmetries acting on bilayer systems, assuming that 
each single layer systems exhibits emergent modular sym- 
metry in its scaling behavior. In the most general case 
a bilayer system has two independent complex conduc- 
tivities, one for each layer, plus an inter-layer conductiv- 
ity. A genus two modular group S'p(4, Z) then acts on a 
complex two-dimensional matrix of conductivities (with 
positive imaginary part). If the bilayer system consists 
of identical layers, then the scaling behavior is expected 
to be the same as in a mono-layer system, and the ex- 
perimentally accessible total conductivity a for the bi- 
layer is just the double of the conductivity a for a sin- 
gle layep^Sl. if the scaling properties of a system with 
a single fermionic degree of freedom is determined by 
Tt, then a = 2(j implies that the scaling of a bilayer 
system is constrained by Tr. One example of this is ob- 
served when a spin polarized system changes to a spin 
degenerate system, summarized by the symmetry trans- 
mutation Tfi ~^ ^ > T t that is induced by decreasing 
the magnetic fielcP^'^. On the other hand, if the scal- 
ing of a single bosonic degree of freedom is determined 
by Ts, then the scaling of the total conductivity a of 
a bilayer system should respect a symmetry generated 
by T^{a) = a + 4: and S'(ct) = -4/5-. This is verified 
by using the duality transformations and S acting 
on a = a/2. Although S is not a modular transforma- 
tion, the substitution a — > cr should produce a modular 
scaling diagram. Ref. [31] compares these predictions with 
the fixed point structure observed in multi-component 
quantum Hall experiments. While the experimental sit- 
uations is still quite murky, there does not appear to be 
any inconsistency. 

DolarPSl has argued that one should expect the dual- 
ity group to be r(2) in systems with a small but non- 
vanishing difference between adjacent Landau levels of 
different flavors (spin, layer, etc.), so that there is mix- 
ing between the Landau levels. The degenerate bilayer 
case is the limit of exact degeneracy, and the single layer 
case the limit where there is no mixing between Lan- 
dau levels. The transformation ST^S is the 27r-shift (or 
flux attachment) in the statistics of the charge carriera^Sl, 
and should therefore always be present in two dimen- 
sions. The translation counts the degeneracy of 
Landau levels, since the Hall conductivity transforms as 
T'^ia^) = aH + N. The symmetry is therefore r(2) for 
fermionic charge carriers if one assumes that the chemical 
potential cannot be tuned in the gap between the differ- 
ent flavors of Landau levels. If both particles and vortices 
carry charge the situation is more complicated, but since 
ST'^S is the periodicity in the statistics angle and one 
expects that there is no mixing between Landau levels of 
a single flavor of charge carriers, one should expect that 



the charge transport phenomena in two dimensions is at 
least r(2)-symmetric. 

There have been two other attempts to construct fam- 
ilies of interpolating /3-functions for the QHE?^"^. They 
have both retained the original idea from Ref.i2j elab- 
orated in Ref.O that the physical (contravariant) /3- 
function should be automorphic under modular trans- 
formations with weight w = —2. This is sufflcient to give 
pictures of automorphic vector fields that look similar to 
ours, since any automorphic meromorphic function will 
by definition have this property. They do not, however, 
consider the physical properties of critical points, nor the 
experimental constraints discussed in the introduction. 
Consequently, the conjectured /3-functions are not well 
motivated, nor do they appear to have any reasonable 
physical interpretation. 

DolarP^EH considered interpolations between the max- 
imal symmetries Tx {X = R, S, T) by postulating a fam- 
ily of meromorphic functions that is conjectured to be 
the /3-function for the RG flow. However, these "scaling 
functions" have no physical foundation or interpretation. 
In particular, they have poles where there should be crit- 
ical points, i.e., where a physical /3-function must vanish. 
Critical exponents are therefore ill defined, and there is 
no physical scaling. 

Georgelin et al.!^ have also considered the symmetry 
r(2), and proposed an interpolating /3-function that in 
our notation is /3 cx (a A — l)/^!. They conclude (cor- 
rectly) that the location of the zero of this function is 
not predicted by r(2) symmetry, but they seem ambiva- 
lent about the order of the transition point, and since the 
function is holomorphic it cannot be physical. They do 
not consider the family of functions under deformations 
of a, and do not discuss the points of enhanced symme- 
tries. 

This is in sharp contrast to the idea proposed in Ref. 
1131 and pursued here. Our /3-function is derived from 
a physically motivated and physically sensible potential, 
and therefore by construction is well behaved everywhere 
on the interior of parameter space. This includes the 
critical points where it has simple zeros, and therefore 
well defined scaling and critical exponents. 



V. SUMMARY 

A phenomenological analysis using global discrete modu- 
lar symmetries of RG flows has previously been success- 
fully applied to the spin-polarized QHE^. Following a 
proposal originally made in Ref.ll3( we have here gener- 
alized this to more complicated situations where spin or 
other "flavors" of charge carriers are relevant. In some 
simple cases, including the spin-polarized system, maxi- 
mal subgroups of the modular group PSL(2, Z) appear to 
account for the data, while other systems are less sym- 
metric. 

Motivated by the C-theorem, which offers an appealing 
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physical interpretation of the RG potential, and insisting 
on some reasonable physical properties of this potential, 
we are led to propose the simplest possible ansatz for 
a family of potentials $a that can interpolate between 
these symmetries (eq. ([t])). $a is parametrized by a sin- 
gle number a, up to an undetermined normalization. The 
potential is always symmetric under the main congruence 
group at level two, and when a takes certain values this 
symmetry is enhanced to one of the maximal subgroups 
of the modular group. The covariant RG /^-function is a 
holomorphic vector field derived from $a, and we have 
compared the geometry of this gradient flow with avail- 
able temperature driven scaling data. 

Until substantial progress is made in deriving low en- 
ergy effective field theories from the micro-physics of 
quantum Hall systems, a remains an uncomputable, non- 
universal, sample specific number that parametrizes our 
ignorance about these models. This is no different from 
other effective field theories, and as usual this is not an 
obstruction to progress since the parameter values can 
always be simply extracted from a few experiments if the 
model is correct. In our case the value of a is most eas- 
ily determined from experiment by finding the location 
of a quantum critical point, i.e., an unstable zero of the 
/3-function given by a saddle point of the RG potential. 
The surprise is that a single parameter so far has suf- 
ficed to account for almost all scaling data, in systems 
with widely different micro-physics but similar emergent 
modular symmetries. 

Almost all the data (Refs. l^Bl and Figs. 6-10) are con- 
sistent with a e K, which together with r(2) symmetry 



implies a generalized semicircle law. The semi-circle law 
has in some cases been confirmed by experiment to be 
remarkably accurat^^. 

We have argued that this family exhausts all possible 
RG fiows automorphic under the modular group, given 
some mild constraints on the fixed point structure (the 
plateaux and quantum critical points) of the global phase 
diagram of the quantum Hall system. There are two 
other modular symmetries at level 2, but neither r(l) 
nor Tp have physical potentials. Fortunately, there is 
also no experimental evidence that these symmetries are 
of interest. 

This is one example of a remarkable convergence of 
quantum Hall physics and modular mathematics. Hall 
quantization is, in fact, an automatic and unavoidable 
consequence of modular symmetry. The geometry of 
holomorphic modular symmetries just barely admits the 
existence of functions and forms that we need in order 
to give a quantitative description of the geometry of ob- 
served RG fiows, but no more, and unphysical holomor- 
phic /3-functions are prohibited by the symmetry. 

The rigidity of holomorphic modular symmetries al- 
lows us to harvest an infinite number of robust predic- 
tions: detailed quantitative predictions that are so con- 
straining that it seems unlikely that they can be satisfied 
by any real physical system. Nevertheless, almost all 
low temperature quantum Hall scaling data are so far in 
agreement with these predictions, as shown in Ref.i4jfor 
the spin-polarized case and for the fully spin-degenerate 
case, and for some multi-component systems here. 
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